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Abstract 

We briefly explain the notion of exclusion statistics and in particular discuss 
the concept of an ideal exclusion statistics gas. We then review a recent 
work where it is demonstrated that a two-dimensional Bose gas with repulsive 
delta function interactions obeys ideal exclusion statistics, with a fractional 
parameter related to the interaction strength. 
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I. INTRODUCTION 



Low- dimensional electron systems have been extensively studied during the last decades, 
both due to their practical importance and because of their exotic quantum mechanical 
properties. Prominent examples are the 2D heteroj unctions that at low temperature, and 
in high magnetic fields, exhibit the quantum Hall effect and carbon nanotubes, which 
are quasi Id and (for certain geometries) are expected to behave as Luttinger liquids [§]. 

At the same time there has been a tremendous advance in the techniques for trapping 
and cooling atoms, which has opened new vistas for the study of quantum gases, and in 
particular it is now possible to create quasi-two dimensional atomic Bose gases. 

A very interesting property of certain low- dimensional systems is that the elementary 
excitations obey exotic quantum statistics. The most celebrated example is the quantum 
Hall effect, where the Laughlin quasi particles are anyons ||, i.e. particles whose quantum 
statistics "interpolates" between that of bosons and fermions in the sense that their many- 
body wave function picks up a general phase exp (i8) (rather than just 1 or —1) under 
particle exchange. 

Some years ago, Haldane proposed another type of generalized quantum statistics, so- 
called (fractional) exclusion statistics (FES) 0]. Rather than modifying the many-body 
exchange phases, he generalized the Pauli principle, i.e. introduced new rules for occupying 
single-particle quantum states. The basic idea of FES is that adding a number of particles, 
AiV, to a system, blocks Ad of the states available for the next particle according to the 
linear relation Ad = —gAN. This corresponds to a repulsion between the particles in 
phase space, and only very special types of interactions give rise to this type of exclusion of 
single particle states. One of Haldane's original examples was a one-dimensional spin chain; 
moreover, he showed that the Laughlin quasi particles (or generally anyons in a strong 
magnetic field, i.e. confined to the lowest Landau level) also obey exclusion statistics with 
an exclusion parameter g = O/n. While anyons for topological reasons can exist only in one 
or two dimensions, there is no such restriction on exclusion statistics particles. Nevertheless, 
until some time ago, all established examples of FES were either strictly one-dimensional 
or, as anyons in the lowest Landau level, effectively one-dimensional |7]-|9| . 

In a recent work Jl0| we established that within a certain (experimentally feasible) pa- 
rameter range a (quasi-) two-dimensional Bose gas with repulsive delta function interaction 
is in fact an ideal FES gas where the exclusion parameter is determined by the interaction 
strength. 

The model Hamiltonian, given by 

N ( v 2 \ nh 2 N 

1=1 \ / i<j 

has previously been used to describe atoms in Bose condensation experiments using highly 
asymmetric traps fll|jl2|l . The particular parametrization of the delta function potential 
is chosen as to reproduce the s-wave scattering phase-shift in three dimensions, and the 
dimensionless coupling g is given by 
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where a is the (three dimensional) scattering length and l z the out of plane extension of the 
asymmetric trap, which is harmonic in the transverse direction with a frequency u = Ti/ml 2 z 
|TT| , |T3|| . We shall assume that the temperature is sufficiently high above the transition 



temperature that the only relevant mean field is the density, n, and that the fluctuations 
are small enough to be ignored. We shall return to these issues in the discussion (section 

Our study was motivated by the previous observation [|T^l|T3| that, in a Thomas- Fermi 



approximation, a two dimensional Fermi or Bose gas with short range repulsive interactions 
has the same energy and number density as an ideal FES gas (treated in the same approxi- 
mation). To reach our conclusion, we analyzed the statistical mechanics of the delta function 
gas, both by using a mean field approximation for the energy levels, and by directly relating 
the thermodynamics to the quantum mechanical scattering phase shifts. 

In order to make this presentation reasonably self contained, section |lj gives a brief 
introduction to exclusion statistics and in particular to ideal FES gases. Section [TTI| contains 
the mean field argument for the equivalence of (|IJ) to an ideal FES gas, and in section [IV 



we 



show that this result is consistent with a full quantum mechanical treatment of the model 
within a certain parameter regime. Finally, we discuss our results and point out some open 
questions. 



II. EXCLUSION STATISTICS 

As already mentioned, exclusion statistics with statistics parameter g is defined by the 
relation M 

Ad(N) = —gAN (3) 

stating that d(N), the dimension of the single-particle Hilbert space when N — 1 states are 
occupied, decreases by g whenever a particle is added to the system. Obviously, g = 1 
corresponds to fermions (Pauli blocking), while g = describes Bose statistics. In general, 
g can be any rational or integer number. The corresponding statistical weight, i.e. the 
number of ways in which a given number of particles can be distributed over a given number 



of states, is [15.16 



[D-{g-l)(N-l)]\ 

N\[D-gN + g-l]\ K J 

where D = d(l). Eq.(|4]) is seen to interpolate between the familiar expressions for bosons 
(g = 0) and fermions (g = 1) in accordance with Eq.(|3|). 

An ideal FES gas is defined by demanding that this generalized exclusion principle 



applies locally in energy space i.e. to a set of degenerate or almost degenerate energy levels. 
In the case of the quantum Hall effect, g is, for topological reasons always rational while 
in our case, being related to the interaction strength, it can take any real value. For non 
rational g, W in (J|) is not integer, and has no combinatorical interpretation, but is rather 
used to define the concept of FES. 

In direct analogy with the case of bosons and fermions, the statistical distribution func- 
tion is derived from (|j) by maximizing the non-equilibrium entropy S = -^Ei;' 11 ^! 
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where the sum is over all energy levels e k} subject to the constraints of keeping the total 
energy and number of particles fixed [jO^U. The entropy is found to be 

S = ^2 D k [(1 - (g - l)n k ) In (1 - (g - l)n k ) - n k \nn k - (1 - gn k ) In (1 - gn k )} , (5) 
k 

where n k = N k /D k , and the resulting statistical distribution function (or occupancy factor) 
can be expressed as [|15HT7| n k = (w k + g)^ 1 where 

w l( l +w k ) 1 ~ g = exp((3(e k -//)), (6) 

with (3 = 1/ksT. The special cases g — 0, 1 correspond to the familiar distribution functions 
of ideal bosons and fermions, respectively. At zero temperature n k is a step function, 

n° K = ~0(e F -e k ), (7) 

where ep is the Fermi energy. This means that in the ground state, there are on average 



1/g particles in each original quantum state. In section |T|, we shall derive the expression 
(§) for the entropy starting from the 2D delta function gas ([!]). 

The thermodynamics of exclusion statistics has been studied in a number of papers 
ji6HT9|. A result that is particularly interesting for our purposes, is the virial expansion 
(i.e. the expansion of the pressure in powers of the density) of an ideal FES gas. As shown 



by Isakov et al. |L7|], the equation of state can be evaluated exactly in the special case D = a 
where D is the number of dimensions and a the power of the dispersion, e(p) oc p a . Obviously 
the case we are interested in, a non-relativistic two dimensional gas with e(p) = p 2 /2m, falls 
into this class. Thus, for D = a = 2, the virial expansion of an ideal FES gas takes the 
form 

P 1 1 00 R 

^ n+ I(2 9 -l)A^g^(4«7 (8) 



where At = y 2nh 2 (3 / m is the thermal wavelength. Here it is also worth mentioning that for 
ideal quantum gases in 2D (bosons, fermions or "exclusons"), the pressure equals the energy 
density S. Note that all statistics dependence in (§) lies in the second virial coefficient, that 
interpolates linearly between the familiar results for ideal bosons and fermions in 2D 
In section |V| we shall derive such a linear relation for the 2D delta function gas. 



III. MEAN FIELD THEORY 

Before going into a statistical mechanical analysis of our model, we give a simple ther- 
modynamic mean field argument for why we expect exclusion statistics. For simplicity we 
consider the case with a constant external potential V, so that the density, n, is also con- 
stant. The crucial observation is then that in a mean field approximation, and for a fixed 
number of particles, the interaction energy term in ([]]) just amounts to a constant shift of 
the energy density, 
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irh 2 2 

£ = £-FreeBose H #™ • (10) 

m 

Expressing this in terms of the thermal wavelength, we see that this shift is identical to 
that in the FES virial expansion Eq.@ if we identify the interaction strength g with the 
exclusion statistics parameter. In other words, our interacting Bose system (|1|) has the same 
(mean field) energy density as an ideal FES gas with statistics g. 

Let us now consider the statistical mechanics of ([I]), and assume that the potential V is 
slowly varying compared with the thermal wavelength, At- We then divide the system into 
cells of area b 2 , where At <C b <C |VV/V|, and study the statistical mechanics in each cell. 
In a mean field approximation, the one-body Hamiltonian in the cell I becomes 

n 2 2nfi 2 

Hi = ^- + V(fi) + gn(n) , (11) 

2m m 

with the corresponding energy e e , where we approximate the potential in the cell with the 
constant V(rt) with fi the position of the center of the cell. Also n{r^) = Ng/b 2 is the mean 
number density in the cell £, with Ng ^> 1 the corresponding average number of particles. 
[] We want to demonstrate that the number of available single-particle states in a box £ 
decreases with the number of particles as in Haldane's definition of FES, Eq.(|3]). To this 
end, we write the quantized momenta in the box as p = (2it%/b)l with l x and l y integers. 
The total number of states di below some energy e e corresponding to a momentum p, is 



di = irl . Using (11) to eliminate p, one finds 



de = ^[e e -V(r e )]-gN e (12) 
= d/,BoBB - gN t . (13) 

Indeed this result indicates exclusion statistics. Note that an important ingredient of the 
derivation of de is that the dimension matches the power of the dispersion, in this case 
D = a = 2. 

Eq. (|T2|) describes the dependence of the total number of states below some energy on the 
total number of particles in the box. But, as explained in section ||, for an ideal FES gas, the 
Haldane exclusion property must hold locally in energy space, so the state counting result 
(12]) is a necessary but not sufficient condition. The strategy to prove that the box Hamil- 
tonian ([□]) indeed describes an ideal FES gas, is to first map the interacting system onto 
a non-interacting one in terms of renormalized quantum numbers (pseudo-energies) which 
display the exclusion statistics property. Then we can use standard counting arguments to 
arrive at the expression (R) for the entropy of an ideal FES gas. 



In going from (|lj) to (11) it is important to correctly incorporate the effect of Bose statistics [21]. 
There is an extra factor of 2 in the interaction term in the Heisenberg equation for the quantum 
field operator as compared to the similar looking Gross-Pitaevskii equation for a classical Bose 
field. This also implies an extra factor of two in the mean field one-body Hamiltonian (]Tl|), and 
here we differ from the treatment in [13|. 
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The starting point is the total energy, Eg which is not simply the sum of the one particle 
energies e\ of (|TT1), but given by 



e i gn{r e ) 

m 



(14) 



where the last term compensates for the double counting of the interaction energy. Re- 
calling from above that the number of states below a given energy is proportional to the 
corresponding kinetic energy, we can rewrite Eq.(|14D as 

Nt ( 2nh 2 a \ 

^ = ^ i { m) + ^ [h+ 2 Ne] ) ' (15) 

where ki = dfj is the number of available single-particle states below the energy of particle 
i, and we choose to order the labelling of particles in a given microstate such that < 
k\ < k 2 ■ ■ ■ < k Nr Although the energy levels (the integers {h}) are not equally spaced, we 
assume that the box is large enough for this effect to be negligible. 
Next introduce the quantities ki by 

k i = k + gY J Kk-~h) , (16) 



with 6{x) being the step function, and rewrite the energy as a sum of "pseudo-energies", e 



1 1 



Et = E~4 (17) 



i=i 



2nh 2 ~ 



Note that, expressed in terms of these pseudo-energies, the system is formally non- 
interacting. The interaction term ~ gN^ has been absorbed in the renormalized quantum 
numbers fcj. The exclusion properties of this system are now manifest since (|T6| ) implies that 
the pseudo-energies must satisfy, 

~4 +1 >4 + ^g. (18) 



As mentioned above, the relation (J12[) holds only because the kinetic energy and the 
number density scale as the same power of the cell size b. This is true for the present 
case of particles with quadratic dispersion in two dimensions, but also for particles in one 
dimension with linear dispersion, which is the case for anyons in the lowest Landau level, 
or equivalently, chiral particles on a circle with an iV 2 type interaction ||. In fact, the two 
models studied in this paper and in can be exactly mapped onto each other by identifying 
the ki in (0) with the "pseudo-momenta" introduced in 0. 

The proof that the exclusion property ([18]) corresponds to an ideal FES gas as defined 
in [|15[], amounts to showing that the (non-equilibrium) entropy of our system equals that of 



an ideal FES gas. This derivation is a straightforward modification of the one given in [22 
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for a multispecies system in the fermionic representation: Going to a continuum description 
with "momenta" pi and "pseudo-momenta" pi defined as 



'2%h\ 2 . {2nh^ 2 



Pi = ki 'i Pi = \ ~b~J ki > ^ 

replacing the sums over k and k by integrals in the usual way and denoting the corresponding 
particle densities in momentum space by v(jp) and p(p), respectively, one finds the continuum 
version of Eq. (|16|) , 

P — P + 9 J dp'p(p')6(p — p) . (20) 
Furthermore one has to demand conservation of the number of particles when changing 



variables from p to p, i.e. v(p)dp = p(p)dp. Combining this with Eq . (|20D gives 



Inserting this into the standard expression for the bosonic non-equilibrium entropy, 



S 



~ k {ih) J d P\v^v-{\ + v)hy[l + v)] , (22) 
exactly reproduces the continuum version of the FES entropy (0) 

S = ' k {^hj J d P{p^P-(i-(g-i)p)^(i-(g-i)p) + (i-gp)^(i-gp)] , (23) 

from which all thermodynamics follows. 

IV. BEYOND MEAN FIELD THEORY: QUANTUM MECHANICS 

So far our analysis was entirely in the context of mean field approximations. It is an 
interesting question whether the full quantum problem of a two dimensional gas with a 
delta function interaction also allows a description in terms of exclusion statistics in some 
range of temperatures. In this section we show that this is indeed the case, by computing 
the correction to the second virial coefficient due to the interaction, and showing that in a 
certain parameter regime, the result is consistent with Eq.(^J). To this end, we first calculate 
the quasi-2D scattering solution with a delta function potential to find the s-wave scattering 
phase shift. The corresponding shift in the pressure of the gas can be found from the Beth- 
Uhlenbeck formula |23| which relates the second virial coefficient to the scattering phase 



shifts. Although the interaction naively does not involve any dimensionful parameter, it is 
known that a pure delta function interaction gives rise to short distance singularities and 
requires a renormalization of the interaction strength which introduces a renormalization 



scale and thus breaks scale invariance |p^| . The two-dimensional s-wave scattering phase 
shift 5q is given by p5| 



7 



COt 5n 



I 2. (pl e /f s 



+ - In 

ng 7T 



h 



(24) 



where p is the relative momentum in the two body scattering process, and g is related to 
the (three dimensional) scattering length by (0). Note that the phase shift does depend on 
the momentum p via the renormalization scale, l e J* which up to a numerical factor equals 
the transverse extent l z of the quasi 2d system. The Beth-Uhlenbeck formula then gives the 
following shift in the pressure due to interactions [21 



AP 



in\ T )' 



2kT 



71 



dp e ™ kT — — 
dp 



(25) 



In general, the integral in (|25|) is a function of I z /\t, but in the relevant parameter range, 

a < l z < A T , (26) 

it is approximately constant, as can be seen by the following argument: First partially 
integrate and note that the phase shift vanishes at p — 0. For small g, i.e. a <C l z we 
consider the 0(g) correction to the ideal Bose gas: 



dAP 
dg 



ran 



■(nXq 



dppe~ 



. d5o 
dg 



(27) 



where 



d5o 
dg 



ng' 



1+ \-— + -ln{pl e / f /h) 
\ ng 7T 



(28) 



Next note that the integral has a high momentum cutoff at p ~ due to the exponential. 
This cutoff suppresses contributions from values of p large enough for the logarithm to give a 
sizable correction, since we have to demand l z in order for the system to be effectively 

two-dimensional. For small p, the contribution from the logarithm is suppressed for small 
enough g, i.e. by demanding exp(l/g) ^> Xt/I z - Neglecting the logarithmic dependence 
in the integral in (|27j), one finds ^ ~ — | . Thus, in the parameter range ( p6| ) there is no 
dependence on the scale l z , and the Gaussian integral can be performed to get the following 
result for the leading order in g correction to the pressure, 



AP = -kT(n\ T ) 2 g 



(29) 



in perfect agreement with ([T0|). We can thus conclude that for temperatures and couplings 



in the range (26) the scattering approach is consistent with the mean field approximation 
used in the previous section. [| It is an open question whether there are any corrections to 
higher virial coefficients. 



'This range might be experimentally relevant in the near future [27|. 



S 



V. DISCUSSION 



Since we have explicitly ignored both the possibility of a quantum condensate and of 
pairing fields other than the density, the results of this paper (and those of Ref. can 
not be used for temperatures below or in the vicinity of the Bose condensation transition 
T c . This is true irrespective of whether this transition is of the Kosterlitz-Thouless type or 



not [ 25 1 . Rather our results should be relevant in a temperature regime where the exclusion 
statistics, due to the repulsive interaction, corresponds to a small correction to the ideal 
Bose gas. Moreover we have seen that one has to demand a <C l z <C Ay in order for the 
scattering approach to be consistent with the mean field results. Since this sets a lower limit 
on the thermal wave length it also implies that the mean field approach can not be trusted 
for very rapidly changing potentials even at high temperatures. 

It is an interesting open question whether the quasi particles above a two-dimensional 
Bose condensate also can be described using exclusion statistics. To answer this question 
one would analyze the corresponding statistical mechanics in a more sophisticated mean 
field approximation that includes effects of phase coherence and pairing mean fields ||21|| . 

Acknowledgement: S.V. thanks Prof. Christophe Salomon for a discussion on possible 
experiments. 
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